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ABSTRACT 

According to the recent theoiy of speciation through loss of 
planktotrophy, the pairs of northeastern Atlantic species of 
caenogastropods would be the result of a cladogenetic event in 
which a planktotrophic ancestor gave rise to two geographically 
separated species with different inodes of laiwal development, 
but retained virtually identical tclcoconch characters. This idea 
was proposed as the working h)pothesis for the origin of tlie 
pair of supposed sibling species Rissoa gaciinii Recluz, 1843 
and R. lia (xVIontcrosato, 1884). The present study shows that 
claims of a close moiphological resemblance between the t\\'o 
species are unjustilu'd as a considerable divergence emerged in 
shell geometric moiphometric analysis and in 16S rRNA mito¬ 
chondrial gene. We conclude that R. giichnii and R. lia cannot 
be regarded as sister or ci')ptic species and should no longer be 
considered a planktotrophic/non-planktotropliic pair. 

Additional keywords: Gastropoda, sibling species, laiwal devel¬ 
opment strateg)^ mitochondrial DNA, geometric moiphometiy 


INTRODUCTION 

Many authors (e.g. jablonski and Lutz, 1983) identified 
hv^o main categories of developmental strategies in 
marine invertebrates: plaiikdrotrophy (P), with lanae 
feeding on plankton, and non-planktotrophy (NP), with 
planktonic lamie feeding only on tbc'ii' yolk supply 
Oecithotrophy), or with direct development. 

According to the so-called “shell-apex rule” (Thorson, 
1950), laiwal development in marine gastropods can 
be inferred from obseiwations of the protoeoneh. A 
mnltispiral protoeoneh and a smaller dimension of the 
initial whorl aeeonnts for a planktotrophic developmen- 
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tal mode, whereas a larger paneispiral protoeoneh is 
directly linked to non-planktotrophy. 

There are several examples of nortehasteni Atlantic 
eaenogitstropod pairs of related taxa that possess identi¬ 
cal teleoconehs and differ exehisively in protoeoneh 
characters. Oliverio (1996) provided a list of 28 of these 
P/NP pairs and indicated that their specific status is still 
a matter of debate. 

Verdnin (1986) maintciined that in species of Rissoa 
(Rissoidae) the oeenrrenee of this phenomenon is of 
remarkable extent and nine pairs of Oliverio s list were 
species of this genus. In his revision of the genus, X^erdnin 
(1976, 1982, 1985, 1986) considered differences in dev el¬ 
opmental mode to definitively separate species, thus 
regarding each member of a P/NP pair as a distinct spe¬ 
cies. His opinion vv^as commonly accepted (c.g. Boiichet, 
1989), although some authors have more recently 
suggested the possibility^ of intraspecific variahiht}' in lar- 
v^al dev^elopment (Waren, 1996; Cadee, 1998) or that this 
v^ariahon indicated incipient speciation (Relifeldt, 1968; 
Panico and Patti, 2()()5). 

Oliverio (1996) suggested a meehanisni of speciation 
explaining the origin of P/NP pairs in the NE Atlantic 
region. This involved the modification of lan^al dev^el- 
opment, with one species abandoning planktotro- 
pliic feeding, thus giving rise to anothca' species. From 
the obsen^ation that in the Mediterranean the non- 
planktotrophic mode is more dominant in the East¬ 
ern basin, Oliv^erio (1996) suggested a paleogeographic 
model to explain the reasons of speciation through the 
loss of planktotrophy During glacial periods sea level 
lows caused isolation bctvvwn the Mediterranean and 
tlie Atlantic, and betwuen the Eastern and Western 
Mediterranean basins. The resulting conditions gav e rise 
to factors thouglit to select against planktotrophic laivae 
(lluctuations in tlie energy^ input, restricted areas, and 
liigher predation pressure; Strathmann, 1978a, b). These 
factors may have caused the shift from planktotrophy to 
lecitotrophy resulting in speciation and origin of a P/NP 
pair. 
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\Vithin this scenario, species forming a P/NP pair can 
be considered sibling species (sensu Knowlton, 1986), 
being both ci')p)tic (i.e., difficult to distinguish using the 
traditional moiphological characters) and sister (i.e., 
sharing tlie same ancestor). 

After the examination of a large amount of museum 
material of Rissoa guerinii Recluz, 1843, and Wissoa Via 
(Monterosato), 1884 (Rissoidae), Verduin (1985) 
highlighted a strong resemblance in shell characters 
between these species and concluded that they “may 
often only be identified reliably by their type of apex” 
(Verduin 1986: 14). Their distribution (Verduin, 1985) 
overlaps in the Western Mediterranean, but only 
R. guehmi, occurs in the Atlantic and is absent from the 
Adriatic and Aegean, where R. Ua is present. 

In this work we pro\ade, based on both diy and live- 
collected material, a critical reinteipretation of the two 
taxa and investigate their phylogenetic relationships. 

MATERIALS AND METHODS 

The source of the biological material examined in this 
study is twofold. Snails from field sampling and empty 
shells belonging to the historical collection of Philippe 
Dautzenberg (housed in the Royal Belgian Institute of 
Natural History (RBINS), Brussels) have been used in 
this research. 

Live-collected Material 

Sampling of living material was performed in the 
infralittoral of Santa Tecla, Sicily (Mediterranean, Ionian 
Sea), where both species commonly occur (Sciideri, pers. 
comm.), at 1—5 m depth in Dec. 2005, Apr., Jun., and 
Nov 2006, About 0.03 m^ of the red alga Pterocladiella 
capillacea (S. G. Gmelin) Santelices and Homniersand 
was collected by SGUBA dmng for each sample. Mate¬ 
rial was immersed in seawater and transferred to the 
laboratory. Each sample was di\4ded into 20 subsamples 
that were individually washed in a tank containing 5 1 of 
50% seawater for not more than five minutes. Osmotic 
shock pro\aded caused the vagile fauna to detach from 
the algal thallus and fall in the bottom of the tank, from 
wRere the specimens were quickly collected and placed 
in seawater. Living material was sorted and taxonomi- 
cally determined under a Wild Makroskop M420 stereo¬ 
scopic microscope. Specimens belonging to Rissoa 
guerinii and R. lia were isolated; their sex was deter¬ 
mined by checking the presence of a penis in the right 
side of pallial cavity. Some were placed separately in 
running seawater at 18°G and provided with fresli 
R capillacea thalli, others were immediately presented 
in 80% ethanol. Reared specimens did not suiAive more 
than 6 weeks; empty shells of dead specimens were 
retained and uschI for obsewation. 

Shells and Head-foot: Adult living specimens of l)oth 
R. guerinii and R. lia were placed in a Petri dish with 
seawater. Sliells w^ere held with forceps and the snails 


attempted to crawl extending their foot completely, 
enabling the head-foot to be observed in detail. Images 
were taken and digitized using a Leica DEG 300 EX 
video camera and Leica Application Suite version 2.4.0 
software, Color drawings were also made to better rep¬ 
resent the color pattern of tlie head-foot. 

Traditional Shell Moi'phoinetry: Standard teleo- 
conch (L, M, W^.i, D^.i) and protoconch parameters 
(d and Do) (Verduin, 1982a) were measured on a sample 
of 100 randomly selected shells (from 53 males and 
47 females) from live-collected material of each species. 
The total number of shell whorls (N) was also counted. 
The following standard teleoconch ratios were calcu¬ 
lated: relative height (L/N), slenderness (L/D^.i), rela¬ 
tive aperture height (M /L) and last-whoii height over 
width ratio (W^.^/ Dn.i). 

A principal component analysis (PGA) and a discrimi¬ 
nant analysis (DA) were performed on the dataset 
obtained combining protoconch parameters and 
teleoconch ratios. Tests and plots w^ere implemented by 
SPSS V.15 statistical software (© SPSS Inc., 2006) and 
by SYSTAT statistic software v. 12 (Wilkinson et al, 
1992). 

Geometric Shell Moiphometry: Fifteen shells 
(8 females and 7 males) at terminal growih were ran¬ 
domly selected from live-collected material of each spe¬ 
cies. Shells were observed using a Leica Z16 APO 
stereoscopic microscope, and colour images were taken 
and digitized using a Leica DEG 300 FX video camera 
and Leica Application Suite version 2.4.0 software. The 
shells were always placed in the same position, wdth 
the coiling axis in vertical position and the aperture on 
the same plane as the objective (Garvajal-Rodriguez et ah, 
2005). Using the software tpsDIG2 v. 2.10 (Rohlf, 2007a), 
19 landmarks (LM) were established (Figure 1). LMl is 
the apex of the shell; LM2, LM4 and LM6 are placed on 
the right border of the profile at the beginning of the 
three last complete whorls, LM15, LM17, and LM19 are 
the corresponding landmarks on the left border of the 
profile, LM3, LM5, LM16, and LM18 mark the interme¬ 
diate position respectively between LM2 and LM4, LM4 
and LM6, LM15 and LM17, LMl7, and LM19 along the 
ciiwature of the whorl; LM8 is at tlie lower suture of the 
last complete whorl and LM7 marks tlie intermediate 
position l)Ctween LM6 and LM8 along the ciiwature of 
the whorl. LM9 is the most external position in the upper 
part of the outer lip; LMIO and LM12 are the most exter¬ 
nal positions respectively in the external right and left 
part of tlie outer lip; LMll is the lowest point at the base; 
LM14 is the most external point in the last whorl at the 
left profile of the shell; LM13 is the profile point between 
LM12 and LM14 (closest to LM7). As described in 
Ganujal-Rodriguez ct al, (2005) the matrix of raw coordi¬ 
nates generated by tpsDIG2 was used in tpsRelw v.L45 
(Rohlf, 2007b) to compute shell size (CS), uniform 
(U1 and U2) and non-uniform (several relative warps, 
RWs) shape components for each specimen. Standard 
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Figure 1. Shell of a live-collected R. lia from S. Tecla show¬ 
ing the placement of the 19 landmarks used for geometric 
moiphometric analysis. Scale bar = 1 mm. 


parametric tests were performed on tlie obtained vari¬ 
ables using the SPSS/PC package v. 15.0. 

Raclular Analysis: The shell of 10 live-collccted spec¬ 
imens of each species was removed and the bodies incu¬ 
bated for 3—4 hours at 50°C in 70% KOH solution. After 
complete tissue dissolution, radulae were isolated, rinsed 
in distilled water, and mounted on SEM stubs. Obseiwa- 
tion and pictures were made using a Jeol JSM-6700F 
scanning electron microscope. The cusps of rachidian 
and lateral teeth were counted and frequency histograms 
were drawn to show differences in cusps arrangement 
between the two species. 

Molecular Systeiiiatics: Foit)^-five ethanol preseivcd 
specimens of l\. giierinii and 34 of R. ha were randomly 
selected from live-collected material. The shell of each 
specimen was broken in a mortar and removed. DNA 
extraction plus amplification, purification and sequencing 
of a 337 bp segment of mitochondrial 16S rRNA gene 
were performed as described in Criscione et ah, 2009. 

Sequences obtained were aligned wdth CodonCode 
Aligner v. 1,6.3 (CodonCode Coiporation, Dedham, 
MA) and the alignment refined by eye. For all samples, 


Table 1. List of studied material in Dautzenherg collection 
(RBINS). 


Acronym 

No. of shells 
contain c'd 

Main original label 

Rga 

About 80 

Rissoa giiciiai Reel, van alhina 


K. ouerinii 

Dautz./Rochebonne/dd. 16.IX.05 

Kgb 

About 60 

Rissoa ^iienni Wed. van bipaiiita Dz 


K. ^iiennii 

and Diin/Rochehonne/dd. 16.IX.05 

Rgc 

About 80 

Ris'soa ^uenni Reel, van conspersa 

R. giicnnii 

Dautz. and Duroiichoux/ 


Rochebonne/dd. 16-IX. 05 

Rlma 

Many R. lia 

R {Apiciilaria) lia, Benoit/t)pique de 
Messinal/Monterosato 2.III. 17 

Rlmb 

5 K. lia 

his/R. (Apiciilaria) lia, Benoit/ 
Messina! Playa 

Rip 

8 R. lia 

Rissoa lia Monts./Paulille/Buccpioy 

Rlt 

many K. lia 

Rissoa lia Benoit/Tanger/Pallaiy 1. 7. 
98 


Table 2. GeiiBank accession numbers for the sequences 
used in the analysis. Accession numbers wdth the prefix GU 
were collected for this study. 


Taxon 

AN 

Analysis 

R. giierinii 

from GUI77882 to 
GUI 77921 

NJ, MP, ML, MJ 

R. lia 

from GUI77922 to 
GUI 77955 

NJ, MP, ML, MJ 

R. siinilis 

GU 177963 

NJ, MP, ML, MJ 

R. lahiosa 

AY676117 

NJ, MP, ML, MJ 

R. pai'va 

AF445343 

NJ, MP, ML, MI 

R. aiiriscalpiiim 

GUI 77879 

MJ 

R variahilis 

GU 177880 

MJ 

R. violacca 

GU 177881 

MJ 


both fonvard and reverse strands were analysed. 
Genbank accession numbers for sequences used in the 
analysis are given in Table 2. 

For the sequences generated. Parsimony and Maxi¬ 
mum likelihood trees were obtained using PAUP* v. 
4.04 (Sw^offord, 2003). The program Modeltest version 
3.06 (Posada and Crandall, 1998) w^as employed to select 
IIKY+I model in Maximum Likelihood. In computing 
trees, the option of 1000 bootstrap replicates was 
selected. A reduced median joining network (MJ) 
(Bandelt et ah, 1999) was obtained with the software 
Network w 4.5 (http.7Avww.fluxus-technology.com/ 
NET\V4500.exe). 

Dautzenberg Collection Material 

The relevant material in tlie Dautzenberg collection wtls 
examined using an Olympus® SZXIO stereoscopic 
microscope. After this suiwey, 15 lots were selected as 
on their historical \'ahie as vouchers (Verduin, 1985) 
and the specific determination of the shells contained 
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was verified according to current taxonomy. Clearly 
misclassified specimens in the lots were temporarily 
removed and not further considered. Several specimens 
of each lot were obseiwed using the stereomicroscope 
and digital pictures were taken wdth an Olympus® 
CAM ED IA C-7070 \\^ digital camera. Table 1 contains 
the list of the lots studied, accompanied by the data on 
the original label, the number of specimens present and 
photographed, and their revised species determination. 
An acronym is given to enable identification of the lot in 
the followang text. 

RESULTS 

LiN'E-COLLIXTK d M ate ri al 

Rissoa guerinii Shell: Shell spindle-shaped; apex 
solid, shaip and glossy; whorls 7, apical 2-3 nearly Hat- 
sided, 3-4 youngest ribbed and more convex, \Uth pen¬ 
ultimate whorl generally bulging out and gbnng shell 
characteristic fusiform shape. Shells of females always 
bigger than those of males with same number of whorls. 

Shell ornamentation composed by axial ribs and inter¬ 
vening furrows, spiral ridges and inteiA/ening grooves, 
and gro\\4h lines. Shells devoid of spiral ornamentation 
rarely found. Axial ribs limited to youngest 3-4 whorls, 
11-14 per whorl, robust and prominent, swelling in mid¬ 
dle of whorl and fading toward shell base; normally 
slightly opisthocline, tending to prosocline on body 
whorl near aperture and close to suture being slightly 
flexuous; each rib equal in breadth to interspaces. 


Extremely weak labral rib, often ha\4ng appearance of 
whitish smudge as broad as two axial ribs, starting from 
end of last axial rib and ending before outer lip. Spiral 
ridges with grooves delicately and closely cancellated as 
a result of intersection with grov\4h lines. Gro\\4h lines 
prosocline and running obliquely over ribs and furrows. 
Aperture as described by Fretter and Graham (1978). 
Two distinct shell colour varieties detected: typical and 
“var. conspersa' (Dautzenberg and Durouchoux, 1914). 
T)pical R. guerinii (Figure 2) with ribless white-browm- 
ish or lilac whorls, youngest showing sinuous browm spi¬ 
ral lines corresponding to furrows of ribbed whorls. 
Ribless whorls with fawn or brown background colour 
on furrows among ribs and on shell base (sometimes 
tinged pale lilac up to base of ribs). A checkerboard 
pattern (or a series of zigzag lines) created by back¬ 
ground color and a whitish color occasionally present. 
Ribs always whitish or pale lilac. Peristome pale \4olet, 
throat with browm-lilac band extending on to columella. 

Rissoa guerinii “var. conspersa" (Figure 3) with entire 
shell covered by a uniform chess board pattern or a 
series of zigzag lines (as described for some t)q)ical 
shells), wdtli predominance of a dark brown colour over 
whitish colour. Ribs often encircled by a white line, nor¬ 
mally interrupted in correspondence of inten/ening fur¬ 
rows. Apex, peristome, and inner part of aperture as 
described for L^^ical pattern. 

Rissoa guerinii Heacl-foot: As for shells, two colour 
t}pes detected (Figure 6, 7), \4z. topical R. guerinii 
and R. guerinii “var. conspersa" (Dautzenberg and 



Figures 2-9. Pietuivs ol sliclls and drawings ul soil body parts ol R. guerinii typical (2, 6), R. guerinii “var. conspersa" (3, 7), R. lia 
“var, casla))ea" (4, 8) and R. lia “var. fulva" (5, 9). Scale bar = 1 nnn. Drawings by D, Scuderi. 
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Duroiiclioux, 1914). Foot always whitish and median 
part of sole stained brown, lighter in tx^ncal R. guerinii 
than in “var. conspersar Snout light brown in /i. guerinii 
t)pical and darker brown in "var. ccnispersa” Margins 
of distal portion of snout and remaining part of head 
yellowish in R. gnerinii and light brown in “var. 
conspersaT Cephalic tentacles whitish, but sometimes 
dark brown in “var. conspersa.’^ A wliitish spot behind 
base of cephalic tentacles always present. 

Rissoa lia Shell: Shell conical, apex solid, obtuse and 
opaque; wdiorls 6, all equal-sided and tumid, apical 2 
whorls ribless, others ribbed or ribless (though last whorl 
often ribless), penultimate generally not bulging out rest. 
Sexual dimoiphism as for R. guerinii. Shell ornamenta¬ 
tion of axial ribs and inten ening furrows, spiral ridges 
and inteiwening grooves, and growth lines. No shells 
devoid of spiral ornamentation found. Axial ril)s (w4ien 
present) always 14 per whorl, strong and prominent, with 
same thickness across wdiorl but fading towards base on 
body w^iorl; normally slightly opisthoclinc, tending to be 
prosocline; each rib slightly narrower than inteiwening 
furrow^. No labral rib or whitish smudge before peristome 
present. Spiral ornamentation as for R. guerinii. Aper¬ 
ture oval or D-shaped, peristome not showing sinuses or 
slight projection of inner lip; edge vciy thin, beveled 
internally turning out to form a though thin llange. Outer 
lip arising below pcripheiy of body wdiorl, (somew4iat 
b(4ow^ level at w4hch ribs end), its cunature initially 
slight or sometimes nearly straight. Columella short, 
peristome everted over a groove, no umbilicus present. 
Tw^o distinct color varieties detected, here named 
castanea and fuha (after Monterosato, 1884). Intermedi¬ 
ate specimens rarely found, R lia “var. eastanea'^ 
(Figure 4) \iolet-brownish to dark browm w4th ribs aKvays 
lighter. First two whorls white but sometimes browm to 
dark brown. Peristome \iolet-brownish. Pi. lia “var./n/tY/” 
(Figure 5) uniformly fawn with ribs aKvays lighter or 
even whitish. First two w4iorls fawm or wliite. Peristome 
as in “var. castanea \ but neck with a browai-lilac band 
running from lower border of labral rib to columella, 

RisHoa lia Head-foot: As for R. gnerinii, two different 
color txpes detected remarkably distinct in colour pat¬ 
tern and corresponding to shell colour variety and thus 
named castanea and fitlva (Figures 8, 9). Foot aKvays 
entirely whitish and cephalic tentacles aKvays yellowash 


with a whitish spot behind their base. In Pi. lia “var. 
casianea^' snout dark and its distal part yellowasli, in 
"var. fulva’^ snout yellowasli, often with a short lirownish 
stripe running along median part. 

Traditional Morphoinetiya Among the 6 principal 
components extracted, PCI (55%) and PC2 (26%) 
explained most of the variance ohsen ed (Tabic 3). Com¬ 
ponent matrix (Table 4) and loading plot (Figure 10) 
illustrate the correlation between principal components 
and shell variables. PCI is a strongly positively correlated 
wdth IVN, L/Dn.j and Wn-i/ Dn-u strongly negatively 
correlated with M/L and wealdy negatively correlated 
with (1 and Do. PC2 is not significantly correlated with 
teleoconch ratios hut it is strongly positively correlated 
with d and Do. 

No clear clusters resulted from plotting PCI value of 
each shell against its respective PC2 value (Figure 11). 
Tlic 95% confidence ellipse of Pi. guerinii contains 
almost only positive values of PCI, but both positive and 
negative values of PC’2 whereas that ol R. lia encircles 
mostly negative values ol both PCI and PC2, 

The eigenvalue of discriminant function behvecn the 
two species w^as 1.307, the canonical correlation 0.753 
and the Wilks’ lamlxla (0.433) w^as higlily significant 
(p<0.001). The structure matrix (Table 5) shows the 
cumulative within-groups correlations between discrimi¬ 
nating shell variables and the standardized canonical dis¬ 
criminant function obtained. Table 6 show^s the results of 
the classification statistics obtained using the x aliies of 
the discriminant function of each indi\idual to predict its 
a posteriori species memhership. 

Geometnc Moiqilionietiy : Table 8 shows the per¬ 
centages and a descriptive statistical summaiy ol the 


Table 4. Component matrix for the first hvo principal 
components. 



Component 


1 

2 

D 

-0.271 

0.864 

D(, 

-0.318 

0.848 

tVN 

0.843 

0.181 

DDn.i 

0.914 

0.134 

M/L 

-0.913 

-0.151 

Wn,i/ Dn-i 

0.860 

0.106 


Table 3. Total variance ( 

explained by single principal components. 






Initial Eigenvalues 



Extraction Sums of S(|uartYl Loadings 

Component 

Total 

% of Miriance 

Cumulative % 

Total 

% of \7iriance 

Cumulative % 

1 

3.294 

-54.901 

54.9 

3.294 

54.901 

-54.9 

2 

1.551 

25.852 

80.8 

1.551 

2-5.852 

80.8 

3 

0.432 

7.193 

87.9 




4 

0.-342 

5.704 

93.6 




5 

0.2-35 

3.909 

97.6 




6 

0.147 

2.442 

100.0 
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Figure 10. Loading plot of the sliell v ariables relativ e to the 
first tvv^o principal components. 


Table 6. Classification summaiy for the discriminant analysis. 


Species 


R. gucrinii 

R. lia 

Total 

Count 

R. gucrinii 

93 

7 

100 


R. lia 

14 

86 

100 

% 

R. fficiinii 

93.0 

7.0 

100.0 


R. lia 

14.0 

86.0 

100.0 


Table 7. Multiple regression model testing allometr)^ for the 
non-uniform shell shape v^ariablcs. 


Multiple regression 

Variables in the model 

r^ F 

Name 

Beta 

0.5 20.3*** 

RW2 

-0.648*** 


R\\'9 

0.282** 


< 0.05, < 0.001. 



PCI 

Figure 11. Scattciplot of the first hvo principal components. 
95% confidence ellipses are drawn for each species. DaslK‘d 
ellipse represents R. gucrinii. 


Table 5. Structure matrix showing the loadings ol each shell 
v'ariahle on the discriminant function. The variables are 
ordered by absolutt' size of correlation vvdthin discriminant 
function. 


Diseriniinant Function 

M/I. 

0.717 

IVN 

-0.699 


-0.672 

Wn-V On- 1 

-0.487 

(1 

0.382 

Do 

0.307 


relative score for CS, the two uniform components and 
the first 10 RWs, explaining more than the 91% of 
the overall v^ariation. Table 7 shows tlie residts of the 
allornetric analysis for shell shape measurements 
conducted by step-wise multiple regression analysis for 
centroid size (as dependent variable) and tw^o uniform 
and 29 non-uniform measurements, as independent vari¬ 
ables. The F test of the regression analysis vv^as highly 
significant (p<0.0()l) and tvv^o of the relativ^e w^aips, 
R\V2 and RW9, contributed significantly to the regres¬ 
sion model on the centroid size. Centroid size, uniform 
components and only the first ten relatwe waips were 
considered in the analysis of variance (ANOVA), per¬ 
formed to ev aluate the significance of differences in size 
and shape v^ariables. Shells of K. guchnii and R. lia dif¬ 
fered significantly in CS (p<0.()01), RWl (p<0.05) and 
RW2 (p<().()5). The cumiilatwe results of the analysis are 
shovvii in Tabk^ 8. The significance lev^el obtained for the 
corrected analysis (ANCOVA) vvdth centroid size as 
covariate was not maintained for the difference in RWl, 
but was not affected for RW2. The eigenv alue of the 
stepwise discriminant function betw^cen the two species, 
calculated for the 29 RWs, vv^as 3.199, the canonical cor¬ 
relation 0.873 and Wilks’ lambda (0.238) vv^as liighly sig¬ 
nificant (p<0.001), indicating a good separation bctvv^ecn 
groups. Twelve shape variables contributed to the 
discriminant function (RWl, RW2, RW4, RW5, RW6, 
RW7, RW9, RWK), RW18, RW23, and RW27). 
The standardized coefficient matrix (Table 9) shows 
die relatwe importance of the independent v ariables in 
determining the standardized canonical discriminant 
function. Using the individual vmlucs of the discriminant 
functions to predict a posteriori species memberships, 
26 (86.7%) individuals out of 30, were assigned to the 
correct speck's, leaving only 4 (13.3%) that were errone¬ 
ously classified. Looking at species statistics (Table 10), 
13 (86.7%) specimens of /I giicrinii were correctly 
assigned to this species and only two (13.3%) were 
assigned to R. lia. The same pc'rcentages of correctly/ 
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Table 8. Descriptive statistical snniniaiy and results of AN(WA and ANCOVA for the main shell size and shape variables behveen 
R. giicrinii and R. lia. **P < ().()5, *** P < 0.001, ns == non significant. 


Measure 


CS 

UI 

U2 

HWT 

R\\^2 

R\V3 

R\V4 

RW5 

R\V6 

R\V7 

PAWS 

RW9 

RWTO 

Variance explained 
... Mean 

1020 

-0.001 

O.OOl 

32.2% 

0.012 

19.4% 
-0.011 

14.5% 

-0.002 

8.2% 

-0.004 

4.2% 

-0.003 

3.6% 

0.003 

.3.2% 

-0,003 

2.5% 

0.000 

2.0% 

0.002 

1.9% 

-0.001 

li. ffwrtnu 

SD 

122 

0.014 

0.010 

0.038 

0.016 

0.025 

0.018 

0.012 

0.010 

0.010 

0.010 

0.009 

0.008 

R. Ha 

Mean 

815 

0.001 

-0.001 

-0.012 

0.011 

0.002 

0.004 

0.003 

-0.003 

0.00.3 

0.000 

-0.002 

0.001 

SD 

128 

0.012 

0.()09 

0.016 

0.026 

0.017 

0.012 

0.010 

0.010 

().()()9 

0.008 

0.006 

0.007 

ANOVA 


20.10*** 

0.11 

0.4.5 

4.79** 

7.97** 

()..32 

2.05 

2.42 "" 

1.96”^ 

2.16"^ 

0.06 "" 

3.43 

0.75 "" 

ANCOVA 



1.06’^^ 

0.46 

3.21 

9.96** 

0.20 

1.00'’" 

1.38 

1.19'^" 

1.57“^ 

0.15'^'^ 

1.79 

2.01 


Table 9. Standardized cotdficient inatiix showing the relative 
importance of the sliape variables. 


Discriminant function 

R\V2 

1.501 

RWl 

-1.276 

RW9 

-1.127 

RW5 

0.978 

R\V7 

0.932 

RW4 

0.911 

RW6 

-0.893 

R\V23 

-0.789 

R\V27 

0.716 

R\V18 

-0.685 


Table 10. Sinnman^ of classification results obtained em¬ 
ploying the discriminant function. 

Species 


R. gnerinii 

Pi. lia 

Total 

Count 

71. giieriuii 

13 

2 

15 


Pi. lia 

2 

13 

15 

% 

R. gnerinii 

86.7 

13.3 

100.0 


R. lia 

86.7 

13.3 

100.0 


erroneously classified specimens of R. lia were ohsen ed. 
In Figure 12 the thin plate spline representation allows 
intei*pretation in geometric terms the positive (character¬ 
istic of R. lia) and negative deviations (characteristic of 
R, giieriuii) \nlues for the most significant non uniform 
shape \"ariahle, 11W2, between tlu' two species. 

Radiilar Analysis: No relevant differences emerged 
between the radulae of R. gnerinii (Figure 13) and 
R. lia (Figure 14) in relation to the shape ol* teeth and 
eusps. Haciiidian wnih one median eusp, two pairs of 
lateral cusps and two pairs of basal eusps. Lateral tooth 
with a median primaiy Ingger cusp with respectively 
2-3 inner cusps and 4-5 outer eusps in R. gueriaii or 
2 inner eusps and 3-5 outer cusps in il. lia. Figure 15 
represents the frequency distril)ntion of tlu^ number 
ol eusps of the lateral tooth obsened in the sample 
studied and shows higher \auiability in 71 gucriiiii lor 
this character. 



Figure 12. Thin plate spline representations for R\V2, show¬ 
ing the deformation of the grid for the average values of 
R. giicrhiii (left) and 71 lia (right). 


Molecular Systematics: The topologies of the P and 
ML trees (Figure 16 and 17) w^ere comparable: two main 
clades, corresponding to 71 giieriuii and R. lia, elearly 
separated from eacli other and from the outgroujDS l)y 
high bootstrap \ allies. 

Two indi\idual sequences of 71 giieriuii formed a 
weakly supported elade closely related to the main one 
of R. giKuiiiii. The MJ network (Figure 18) showed two 
consistently distinct (18 differences) clusters of haplo- 
t)pes (corresponding to 71 giieriuii and to R. lia). R. 
giieriuii displayed higher genetic structure than R. lia, 
wdth one haplotype occurring in 29 specimens, 2 occur¬ 
ring in 2 specimens and 11 unique haplotvpes. No single 
haplot)q:)e of 71 lia appeared to be dominant: tlie most 
common ones oeeurred in 11, 8, 5, 5 and 2 specimens; 3 
were unique. Both clusters w^ere consistently distant 
from 5 out ol 6 outgroups, with 71 lia closer than 
R. giieriuii. Unexpectc'dly 71. aiiri.sealpiiiiu revealed veiy 
little difference (1-2) from some R. giieriuii haploLpes. 
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Figures 13-14. SEM photographs or the radular ribbon ol 
R. gticrinii (13) and R lia (14). Scale bar = 10 pin. 


Dautzenberg Collection Material 

Shell Visual Obsei^alioiis: Rissoa guerinii: The shells 
of the lot Rga were eonsiderably worn, younger whorls in 
most uniformly whitish, in some uniformly yellowish; 
older 2-4 whorls always pale violet, Most of the shells of 
Rgb showed worn whitisli/yellowdsh lowermost whorls 
and remaining whorls with the pigmentation described 
for R, guerinii t)pical. Shells of Rgc were lietter pre- 
semnl and showed the chesslioaixl (or zigzag) pattern 
described above for R. guerinii “var. ronspersa\ 

Rissoa Ua: Sliclls ol' Rlma and Rlmb were' topotypes, 
the lormer being also from Monterosato collection. 
Although these shells were rather worn, the two colour 
varieties, descrilied al)Ovc‘ for 7i. lur were still detectable 
on tlie younger whorls but had completely disappeared 
from the whitish oldest whorls, 

Sliells ol Rip and Rlt were worn but a different pat¬ 
tern, with a wliitish spiral band running in the middle ol 
the youngest whorls on a dark brown background. 


12 
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T 

no. of inner cusps 
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1 ^-T 

no. of outer cusps 
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□ R. Ua 


n 



3 3 4 5 


Figure 15. Frequency histograms of distribution of inner 
and outer cusps in the lateral tooth of H guerinii and R lia. 


DISCUSSION 

Shell and Head-foot; While in the early malacologi- 
cal literature satisfaetoiy descriptions of the shell of 
R. guerinii are available (Jeffreys, 1869; Fretter and 
Graham, 1978), descriptive data on R. lia are limited to 
its essential original description (Monterosato, 1884; 
page 139). 

Verduin (1985) provided a reinteiq)retation of the two 
tcLxa, stressing an extremely close morphological resem¬ 
blance beUveen them. However he liased his diagnosis 
only on museum diy material, often in a poor state 
of presemation. As an example, tlie character “punctate 
spiral striae” (Verduin, 1985, pages 112 and 114), 
reported for the shells of both species, should be consid¬ 
ered the result of the deterioration of the subtle reticu¬ 
late pattern, formed by growth lines and spiral ridges and 
easily visible in fresh shells. 

The redescriptions provided here are based on fresh 
shells at terminal growth, which maintain the peculiar 
characters of ornamentation and pigmentation of each 
species, but at the same time are readily matched with 
original museum material. 

Daiitzenbc'i'g and Durouchoux (1914), on the base of 
the material of the lots here named Rga, Rgb, and Rge, 
described four colour v/arieties for R, guerinii: t)pical 
(vviiite with the intemils among ribs brovvm), eonspersa 
(brovvTi background with a chessboard pattern formed by 
veiy small spots), alhina (totally wiiite) and bipaiiita 
(with the first fwe or four whorls dark violet and the 
remaining entirely wiiite). Based on the same material, 
Verduin (1985) ev en suggested the presence of a north¬ 
ern R. guerinii subspecies (showing the latter two color 
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16 


\ 

; 



64 


0 R. smiths 

100 
/I 81 



R. lahiosa 


R. paira 


R. gucnnii 


17 


'' ^ R. smiths 

R.lia 99 

91 



R. lahiosa 


59 



R. giiermii 


R. pan’a 


100 


Figures 16-17. Parsimony (16) and Maximum Likelihood (17) Trees drawn after the analysis ol 16S se{|iu‘nees ol R. giicri}ui and 
R Jia. Cirele radius is proportional to stajiienees fref|neiiey Bootstrap values are given at the left side ol each branch. 


R. violticea 


0 R. giier 'mn (45) 
O /?. /w(34) 

^ outgroups 
^ niedinn vector 




Figure IS. Median Joining Network dra\Mi on the basis of 16S se(|iienees. Number of setpienees employed in brackets. Circles 
radius proportional to haploUpe fretpiency. Dashed lines represent mutational distance higher than 4 (values reported nearby). 


varieties) and a southern one (showing the former ones). 
We eannot see any justification for these claims as the 
varieties alliiiia and hipai-tita are the result of the deteri¬ 
oration of the t)^3ieal pigmentation, as shown by the 
worn shells forming the lots Rga and Rgb. 

/Vlthoiigh his inteipretation of R. ha (based on the obser¬ 
vation of Rhna and Rlmb lots) is generally correct, Verduin 
(1985: 114) added to its species chagnosis the ol)senyition 
that “in a sample from Trapani, Sicily (fig. 25) many colours 


and colour patterns can be diseenied, among which the 
colour pattern cotispcr.sa known in R. gucnnii" 

This obsenntion apparemtK^ suggested to \^erduin 
(1985) the idea to move R. lia from the subgenus 
Apinilarki Monterosato, 1884, to Gouiostonia Villa, 1884 
(the same subgenus of R. ffterinii). We regard his obser¬ 
vation as incorrect and in disagreement watli the original 
description of Monterosato (1884) and watli onr obseiwa- 
tions on the material in this study. We found no e\adence 
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tliat the color pattern conspcrsa belongs to R. lia. 
Although we could not examine the .sample mentioned by 
\^erduin (1985), an examination of the figures he provided 
(25a-h), allow^ed us to assess that these figures are com¬ 
posed by a mixture of shells of 21. lia and R. giieri}ili shells, 
associated because of their similar apex dimensions. 

However, the statement that some R. lia show “the 
presence of broken brownish colour lines which encircle 
the shells” (Verduin, 1985: 114), is in agreement with 
obsenntions made on lots Rip and Rlt and with some 
recent peculiar records of Rissoa sp. from Sardinia 
(Tyrrhenian Sea), showing this same peculiar coloration 
(Fasulo, pers. comm.). Although this coloration is not 
t)pical of the t}pe material of R. lia^ further obseiwation 
on live-collected material from these localities is needed 
to properly address the issue. 

Rissoa species tend to maintain the head-foot colour 
pattern (i.e., position and dimensions of elements as spots 
or blotches), although they may slunv intraspecific varia¬ 
tion in colour intensity (Fretter and Graham, 1978). This 
phenomenon has been employed to support species iden¬ 
tity (R. pawa and R. iiitcrrupta, Waren, 1996; R. fftcrUiii 
and R. pa}ilion)ic}isis, Criscione et al., 2009) or to h)poth- 
esize recent speciation (R. meinbranacea t)pe A and B, 
Relifeldt, 1968; R. aunscalputm type a and b, Colognola 
et ah, 1986). This rule is in agreement wdth the slight 
vnriation in color intensity obseiwed in this study for the 
two varieties of H. gucri/di, but it does not apply to R. lia, 
whose color varieties showed two distinct patterns. 

Traditional Morphoinetiy: Although the present 
moiphometric analysis is restricted to populations from 
a single locality, its results are not only helpful in under¬ 
standing intrapopulation variabilitv' in shell moiphologv; 
but also in drawdng more general conclusions in terms 
of specific differentiation. The differences obsened can 
be expressed in terms of teleoconch and protoconch 
variation, wath the first element being by far the most 
important. Populations of both species show a wide 
variability in protoconch dimensions, which is larger for 
R. giicn}ui. Piissoa lia shows twice the intrapopulation 
variation in teleoconch characters than that observed 
lor R. gtierinii. 

The values of sliell parameters obtained for the two 
species showed a considerable overlap, which is mainly 
dependent on the strong similarity in [)rotoconcii dimen¬ 
sions and not, surprisingly, by the resemblance of the 
teleoconchs. Most shells of 21. giierinii appear to have a 
larger relative size, a more ekmgated aperture and to be 
more sicaidcr than most 21. lia shells. 

discriminant analysis showed that 21. gucrinii and 
21. ha populations can be distinguished mainly for the 
more elongated aperture of 21. gucrinii shells, for their 
larger relative size, and the higher slenderness, The 
importance of })rotoconch is only marginal, 

In summaiy, our results show that teleoconch mor- 
]ilioni(*tiy can be eflicicMitly used to discriminate' between 
R. gucrinii and R. lia, whereas protoconch dimensions are 
less rejn'cscntitivc' of the ovc*rall interspecilic variation 


and less reliable, due to the large overlap shown. This idea 
is in contrast with the vdevv of tlic discriminating power of 
protoconch dimensions (Verduin, 1986). 

Some Methodological Considerations: Cadee 
(1998: 91) critically revised the methodolog)^ used by 
Verduin (1976, 1982b, 1985, 1986) highlighting that the 
measurements of d and Do “depend on the accurate 
vertical position of the shell and the somewhat arbitraiy 
location of the line along which d and Do are measured.” 
Cadee (1998) dealt with the former issue by performing 
repeated measurements of d and Do on specimens of R. 
inenihranacca, whose vertical placement was reiterated 
ten times, and obtained a standard deviation of about 
0.01 mm for both the variables. This value is roughly the 
same as that which separates the mean values of d and 
Do of R. lia and R. guerinii in this study (not shown). 
This means that some of the differences in protoconch 
dimensions might be indeed the consequence of the 
inaccurate measurements, rather than representing real 
variation. The alternative solution of performing mea¬ 
surements on digital photographs has been also recently 
employed (Cri.scione et al., 2009). The second issue is 
that real Rissoa protoconchs often do not correspond to 
the ideal shell apex figured by Verduin (1977; Figure 1), 
but most of them are rather irregular. In protoconchs 
like these there is a large range of possibilities to locate 
the line along wliich d and Do should be measured and 
the choice is largely arbitraiy^ 

Geometnc Moiyjhonietiy : The ANOVA revealed a 
significant larger size (CS) of the shells of R. gucrinii 
compared to those of R. lia, confirming the results 
of visual obsenations and traditional moiphometiy 
ANOVA showed a significant difference (p<0.()5) 
between the two groups in the first (RWl) and the sec¬ 
ond (RW2) non-uniform sliape v^ariables. The signifi¬ 
cance for R\V2 remained unaltered when correcting the 
analysis for CS (ANCOVA), while that of RWl was not 
significant, indicating that the shape difference ex¬ 
plained by that variable w^as dependent on size. This 
means that the two groups differ exclusively on the sec¬ 
ond non-uniform shape variable (RW2), independently 
from the correlation between shape and size, 

Tlie discriminant function calculated from all the 
non-uniform shape variables, was successful in nioipho- 
metiically discriminating the tvvy) groups. R\V^2 vvyis the 
most important variable in determining the distinction 
l)etw^een 21. gucrinii and 21. lia. The mean vnliies ol this 
variable for each of the two groups (positive for 21. lia 
and negative for 21. guerinii) were plotted in a tps repre¬ 
sentation (Figure 12). The plot showed that that variable 
RW2 is a rellection of the most obvious discernable 
shell shape difference. This comprised the more fusi¬ 
form shape of R. giicriiiii, contributed by consistently 
narrower wliorls than those of 21. lia (represented by the 
contraction of the corresponding zone of the grid) and a 
substantially ecpial penultimate whorl. 
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Tlie grid deformation in the plot of R. lia increase 
progressively from the top to the middle-lower as 
expeeted for a more eonical sliape. In the same plot, 
larger deformations of the upper part and the lower of 
the grid for R. ha, account respectively for tlie more 
obtuse apex and the larger aperture oi R. ha. 

Radular Analyses: Although the taxononiical value of 
radular morphology in rissoids is exclusively limited to 
generic level (Ponder, 1985), eomparative analysis of pat¬ 
terns of similarity may be helpful in resolving the com¬ 
plexity of a highly diverse genus siicii as Ris.soa. Along 
with odier charaeters, differences in the size of the cusps 
of marginal tooth luwe been employed to show differ- 
enees betw^een two moiphs of /i. ainiscalpiinn (Colognola 
et ah, 1986) and radular identity has beam uscxl to support 
findings of recent speciation in R. inonhrajiacea morphs 
(Rehfeldt, 1968). 

How^ev^er, detailed reports on the taxonomie value of 
eusps in rissoids are laeking and this makes it difficult to 
give tlie appropriate w^eight to the difference betw^ecai 
R. giierinii and R. ha in the relative number of cusps of 
the lateral tooth. But the pattern emerged appears to be 
rather eonstant and its eventual v^alue as distinctiv^e char¬ 
acter of phylogenetic significance cannot lie excluded. 

Molecular Systeiiiatics: Analysis of mitochondrial 
secjiienees of gastropod taxa, invoked in a recent specia¬ 
tion ev^ent (or an ongoing speciation process), often show 
traces of introgression (e.g., Kirby et al., 1997; Kojima 
et al., 2001). Introgression shovvm by 16S sec|uenees has 
been used to support ongoing speciation in tlie P/NP 
pair of syinpatrie siblings Rissoa aiiriscalpiiun/R. 
itahensis (Panieo and Patti, 2005). In this study, no evi- 
denee of introgression vv^as detc^cted for s)inpatric 
populations of R giicrinii and R. lia, suggesting that 
claims of recent speciation are not acceptable. Ilow^ev^er, 
the evolutionaiy rate of 16S rRNA is not knovvai in 
rissoids and the hypothesis of an earlier cladogenetie 
event, giving rise to R, gnerinii and R. ha, cannot be 
exeludecL The results of our molecmlar analysis also c‘hal- 
lenge the likelihood of this scenario. 

The ML analysis (Figure 17) failed in resolving the 
polytomy resulted for sefpiences of R. giiehuii, R lia 
and one of the outgroups, R. .siaiih.s, as expected if our 
ingroups vv^ere sister species. Before Verdiiin (1985) 
mov ed R. ha to the siibgenus Goniostoaia, a closer rela¬ 
tionship betw^een this species and R. sirnihs vv/as com¬ 
monly accepted (as both members of the subgenus 
Apiciilaria). Our results may be considered in agreement 
with this idea. 

The hvpothesis of the shared origin of/L giierinii and 
R. lia vv'ould reipiire much lower mutational distance 
betw^een the two sister species than betw/een each sister 
species and the outgroups, Our MJ network (Figure 18) 
shows that this is not the ease. R. ha is separattnl from 
the main haplotvpe of R. giierinii by a distance of the 
same order of magnitude than the distance separating 
this species from the other outgroups. 


The significant difference in liaplot)pe diversity of R. 
giierinii and R. ha (shown by the same analysis) may be 
the result of stochastic ev ents in the ev olutionaiy histoiy 
of tliese lineages or may b(‘ alternativ ely related to the 
supposed alternative strategy of lanal dispersal (P/NP). 
Rissoa giierinii displays a higher genetic structure which 
may reflect a higher genetic flovv^, due to its higher 
dispersal capability. Conversely, R. ha shows a lower 
structure vvTich vvx)iild reflect a lower Icwel of gene 
flovw Despite records of R. giierinii planktotrophie 
veligers have been long report(Hl (Lebonr, 1934; Tliiriot- 
Quievreux and Babio, 1975), no evidence for the 
non-planktotrophic strateg)' oi' R. ha are known and sue- 
eessful laboratoiy experiments on spawning are lacking. 
The actual strategv^ of dispersal of the tvvn species is tliiis 
still to be confirmed. 

Suqvrisingly, an outgroup s(H|ueiice of R. aiirisealpiinn 
showed only two differences from tlie main R. giierinii 
haplotype. This short distance fits in the normal intraspe- 
cifie (and even intrapopulation al) variability of R. giierinii 
and represent a clue of previously unsuspected genetic 
similarity between tlie two species. Further investigation 
is required to properly address the issue. 

CONCLUSIONS 

The taxonomic revisions of Verduin (1976, 1977, 1982b, 
1985, 1986), although based exchisivnly on empty shells, 
have representcnl landmarks in the a-taxonomy of 
Rissoa. Despite some evidence of the taxonomic 
unreliability of shell characters in Rissoa (c.g.Wigham, 
1975; Waren, 1996), only a few attempts liav e been made 
to dispute his conclusions. These include criticism of his 
methods (Cadee, 1998) and the utilization of non shell- 
based approaelies (Colognola ef al, 1986; Panieo and 
Patti, 2()()5; Criseione et al., 2009). This trtaid has been 
followed in our vva)rk, vviiich demonstrates the value of 
a synergy between modern techniques applied to the 
analysis of shell charaeters and molecular methods. 
The present study has sliowai that tiu* veiy similar sliell 
moq)holog)^ of K. giierinii and R. lia liides a clear distinc¬ 
tion in other characters, indicating that they should not 
be considered as a eqptic pair ol sister taxa, 
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